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We construct dark solitons in the recently introduced model of the nonlinear dual-core coupler with 
the mutually balanced gain and loss applied to the two cores, which is a realization of parity-time 
symmetry in nonlinear optics. The main issue is stability of the dark solitons. The modulational 
stability of the CW (continuous-wave) background, which supports the dark solitons, is studied 
analytically, and the full stability is investigated in a numerical form, via computation of eigenvalues 
for modes of small perturbations. Stability regions are thus identified in the parameter space of 
the system, and verified in direct simulations. Collisions between stable dark solitons are briefly 
considered too. 



PACS numbers: 42.65.Tg, 11.30.Er 



I. INTRODUCTION 

The concept of the parity-time VT symmetry was orig- 
inally elaborated in the field theory Q], as a general- 
ization of the canonical conservative systems, which are 
based on Hermitian Hamiltonians, for a special case of 
dissipative systems which include exactly balanced and 
spatially separated linear gain and loss. Such systems are 
described by non-Hermitian Hamiltonians, whose Hermi- 
tian and anti-Hermitian parts are spatially even and odd, 
respectively. A distinctive feature of the non-Hermitian 
Hamiltonians, which are subject to the condition of the 
PT symmetry, is the fact that, up to a certain criti- 
cal value of the strength of their anti-Hermitian (dis- 
sipative) part, the spectrum of such Hamiltonians may 
remain purely real (physical). When this occurs a PT- 
symmetric non-Hermitian Hamiltonian can be eventually 
transformed into Hermitian ones by means of similarity 
transformations @. 

In terms of the quantum theory, PT-symmetric sys- 
tems are the settings of theoretical interest. For the re- 
alization of the PT symmetry in real settings, one can 
make use of the fact that the linear propagation equation 
derived for optical beams in the paraxial approximation 
has essentially the same form as the Schrodinger equa- 
tion in quantum mechanics, in one- and two-dimensional 
(ID and 2D) cases alike. In other words, the evolution 
of the wave function of a quantum particle may be emu- 
lated by the transmission of an optical beam, as in both 
cases the wave propagation follows the same principles. 
This fact makes it possible to simulate many quantum- 
mechanical phenomena by means of relatively simple set- 
tings which can be realized in classical optics Q . In this 
vein, the realization of PT-symmetric settings in optical 
systems, which combine spatially symmetric refractive- 
index landscapes and mutually balanced spatially sep- 
arated gain and loss, was proposed in [4j (see also [5[ 
for subsequent early development of optical applications) 
and experimentally demonstrated in [6|. 



Typically, the models amount to the ID or 2D linear 
Schrodinger equations with a complex potential, whose 
real and imaginary parts are, respectively, spatially even 
and odd. Another possibility of the realization of the 
PT-symmetric settings in optics, in the form a dual-core 
coupler, with the mutually balanced gain and loss applied 
to the two cores, was recently proposed in the works I S LZ f 
for stationary regime of light propagation and in |8l4l0l] 
for the bright optical solitons which exist when the arms 
of the coupler obey Kerr nonlinearity. In this last setting, 
the solitons are available in the exact analytical form, 
and their stability boundary can be found analytically 
too @, [lOj] . 

A natural extension of the analysis of the nonlinear 
PT-symmetric systems is to search for stable dark soli- 
tons in them, which is subject of the present work. We 
notice that the dark solitons in a parabolic potential 
with a PT-symmetric non-Hermitian part, where they 
can be considered as the nonlinear modes, bifurcating 
from the first excited state of the linear PT-symmetric 
parabolic pot ential, have recently been addressed in the 
literature (llj . 

An alternative natural setting for the consideration of 
dark solitons in 'PT-symmetric optical systems is pro- 
vided by the above-mentioned dual-core system. As well 
as a broad class of other solutions, dark solitons in this 
system can be easily found in an exact form Q, the ac- 
tual problem being the analysis of their stability and 
interactions. The model is introduced in Sec. II, and 
the modulational stability of the CW (continuous-wave) 
background, supporting the dark solitons, which is a nec- 
essary condition for their stability, is investigated in an 
analytical form in Sec. III. The mathematical framework 
for the full analysis of the dark-soliton stability is in- 
troduced in Sec. IV, and numerical results, which can 
be summarized in the form of stability diagrams for the 
PT-symmetric dark solitons, are reported in Section V. 
Collisions between dark solitons are briefly considered in 
Sec. V too. The paper is concluded by Sec. VI. 
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II. THE MODEL 

We start with the system of equations for scaled field 
variables q\X- 

.dqi d 2 qi i . .n I 1 2 \ ■ /i \ 

i--q^ = ~-q^i + \xiw\ + xm\ ) qi + *7i«i - 12, (ia) 

.dq 2 d 2 q 2 , . , 2 , ,2\ ■ s 
*— = — + (xl<7l| + X2M J 92 - «72<?2 - gi,(lb) 



Here the linear-coupling constant is scaled to be one, pos- 
itive coefficients 71 and 72 account for the gain and loss, 
respectively, in the two cores, while \ an d Xi,2 are real 
coefficients of cross-phase modulation (XPM) and self- 
phase modulation (SPM). 

Since the subject of the work is the existence and dy- 
namics of dark solitons, it is first necessary to address the 
existence and modulational stability of the carrier-wave 
(CW) background, i.e., solutions in the form of 



qi,2 (z, x) = uij exp (—ibz) . 



(2) 



with complex amplitudes 1/4.2 and real propagation con- 
stant b. The substitution of this into Eqs. (fTJ) yields 



i„ 12 _ I71-72IV1-7172 / 7i^7 9 

' Jl |72(X1-X)+71(X-X2)IV 1 3 ' 3 ~ ' ' 

(3) 

while the relative phase, 6 = argu 2 — arg it! , is deter- 
mined by relation 



tan* = [20(7i- 7a ) -1] 



__\/7l72_ 

y/1 - 7172 ' 



(4) 



where 0(x) is the Heaviside's step function. The propa- 
gation constant of this solution is 



b = 



cos S 



7iX2 



tIxi 



V7172 72 (xi - x) + 7i (x - X2) ' 



(5) 



Note that, according to Eq. @, the CW amplitudes in 
the two components are related by |zi2| 2 /|ui| 2 = 71/72, 
which implies the balance between the gain and loss 
in the CW state. Further, it follows from Eq. (O 
that the background amplitudes have a singularity at 
72/71 = (X ~ X2)/(x ~ Xi) 7^ 1, and this solution ex- 
ists only at < 7172 < 1. This last condition has simple 
physical explanation: it requires the gain (dissipation) in 
an arm to be small enough for being compensated by the 
energy flow from the other arm with dissipation (gain), 
the flow being limited by the strength of linear coupling 
(responsible for the power transfer between the arms) 
which in our case is normalized to one. 

In what follows we concentrate on the case of the VT- 
symmetry, with 71 = 72 = 7- Then, it follows from 
Eq. ([3]) that the nonzero CW background may exist only 
with symmetric SPM coefficients, \i = X2, and for 7 < 1, 
hence it is convenient to define 7 = sin<5, with < <5 < tt, 



and rewrite Eqs. (fTJ) as 
.dqi d 2 qi 



+ (Xikil 2 +xM 2 )<7i +ism(5)q 1 - q 2 , 

(6a) 

l-j^ = ~~q^2 + (xM 2 + XiM 2 )<72-isin(% 2 - qi- 

(6b) 



dz dx 2 
.dqi d 2 q 2 



In comparison with the model of the 'PT-symmetric 
dual-core fiber, which was introduced in (g, |9(, Eqs. (j6]) 
include the XPM terms, which implies a non-negligible 
overlap between transverse modes supported by the two 
cores. Recently, it was demonstrated that, in comparison 
with the well-known results for the SPM-nonlinear dual- 
core system with the purely linear coupling [l2| , the addi- 
tion of the XPM terms essentially affects th e sy mmetry- 
breaking transformations of bright solitons [131 ] and pat- 
terns in the form of domain walls [14] in the conservative 
nonlinear coupler, whose model amounts to Eqs. © with 
5 = 0. 



III. MODULATIONAL STABILITY OF THE CW 
BACKGROUND 

CW solutions of Eqs. © with equal amplitudes follow 
from expressions (|3|): 

qf } = pexp [i{-iy(S/2) - ibz] , b = p 2 { X i+x)-cosS. 

(7) 

Here j — 1,2, and components have phase mismatch 6 
imposed by the gain- loss coefficient. 

To analyze the modulational stability of the CW ([?]). 
we use the standard ansatz with arbitrary real perturba- 
tion wavenumber k, the corresponding eigenvalue, /3, and 
infinitesimal perturbation amplitudes, rjj, vf 



1] = P 



j(-tys/2 



-ibz 



(8) 

Then, two branches f3 = (3i, 2 (k) of the linear excitations 
are readily found as 



fa(k) = ±fcxA 2 + V(xi+x), 



(9) 



h{k) ee ±y/\k 2 + 2 cos 6} [k 2 + 2 cos ,5 + 2p 2 ( X i - *)]. 

(10) 

From ([9]) it follows that, for the stability of the back- 
ground, one has to require 



Xi+X>0, 



(11) 



the constraint which is also necessary for the modula- 
tional stability of the CW background in the conservative 
system (<5 = 0), and which is imposed in what follows. 
Equation ([9]) gives rise to two other conditions for the 
modulational stability, 



cosc$>0, i.e., < S < ir/2; 



(12) 
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(Xi - x)p 2 + cosS > 0. 



(13) IV. STATIONARY DARK SOLITONS AND 

THEIR LINEAR STABILITY 
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FIG. 1: Domains of the modulational stability (dashed) and 
instability (white) in the (6, xi) plane for fixed 7 = sin S = 
0.7 (a), and in the (6,7) plane for fixed xi — 0-5 (b), or 
Xi = —0.3 (c). In all the panels x — !• in panel (a) the left 
edge corresponds to the limit form xi + X = °f condition 
(jllj l. while the bottom edge is given by b = — cos(<5). The 
inset in panel (a) presents the stability domain in the (p 2 , xi) 
plane. 



In the case of self-focusing SPM, x < 0, the stability 
domain is determined by Eqs. (fTTj) and (|T2|) [if these two 
conditions are met, Eq. (|13[) is satisfied automatically]. 
This, in particular, means that the stability of the CW 
does not depend on its amplitude p, being determined 
solely by the interplay between the SPM and XPM coef- 
ficients. 

The situation is qualitatively different for the defocus- 
ing SPM, x > 0. Now, one can distinguish the two dis- 
tinct cases. First, if xi > X [this domain is located to the 
right from the vertical dashed line in Fig. [Ha)], then Eq. 
([TBI is reduced to Eq. ([l~2"j) , thus giving nothing new, the 
background being stable at any amplitude p 2 . If, how- 
ever, xi < X [ m Fig. [TJa), this is the domain to the left 
from the vertical dashed line], then, for the stability of 
the background one needs p 2 < p^ax = C0S( V(x — Xi) or > 
equivalently, b < 2\i cos 6/ (x — Xi)- 

In this situation (i.e., when xi < x) the increase of the 
gain-loss coefficient, (i.e. of S), results in narrowing the 
modulational stability domain, which collapses at 6 = 
ir/2 (7 = 1) as shown in Figs. [JJb,c). The limit case of 
Xi = — X deserves special consideration, since in this case 
system (JSJ becomes effectively linear for equal amplitudes 
of components \q\\ = \q2\ = p. As a result, propagation 
constant b — — cos S does not depend on p. Here CW is 
stable when p 2 < p* 2 = (cos<5)/(2x). In the same time, 
value p* defines a global stability threshold: if p < p*, 
the CW background is stable at any value of xi > — X 
[see the inset in Fig.QJa)]. 



Turning to the study of the dark-soliton solutions, we 
focus on the situation when both components have the 
same intensity profile, i.e., 

q J (x,z)=u(x,z)e«- i y s / 2 (3 = 1,2), (14) 

and thus reduce Eqs. © to the standard nonlinear 
Schrodinger equation, 

i^ = --^ + (Xi+x)\u\ 2 u-cos{5)u, (15) 
whose dark soliton solution is commonly known [l5| : 

^(Xi+X) 

Here b is given by Eq. ([7]), and real parameters v and 
w, which determine the "velocity" (in fact, the spatial 
tilt) and the depth of the soliton, are linked by relation 
w 2 + v 2 = 2(xi + x)p 2 - 

Below we focus on the fundamental dark soliton with 
zero velocity v = (alias the black soliton), u s {x,z) = 
uo(x)e~ lbz where 



uq(x) — ptanh I pJ ^ 1 — ^ x 



(17) 



To address its stability, we first notice that the CW 
background must be modulationally stable, hence the 
parameters to be considered are limited by constraints 
(fTT|) -(fT3 |) . Further, to study the linear stability of the 
entire dark soliton (|17p , we adopt the perturbation solu- 
tion as 

qj (x, z) = [u Q (x) + ufa, z) + iu'fa, z)] e *(-i) 3 V2-«- 

(18) 

with infinitesimal perturbation amplitudes u[ 2 (x, z) and 
u'{ 2 (x,z)\. Then, substituting expressions (|18l) into 
Eq. (|15p . we end up with the eigenvalue problem: 

(9 11 

— = Cu, u = col {u'^u", U 2 ,U2}, (19) 
oz 



with operators 



C = 



sin 8 


L_ 


— sin 5 


— cos 8 


-L+ 


sin 6 


-L 


— sin S 


sin 5 


— cos S 


— sin 5 


L_ 


-L 


sin 5 


-L+ 


— sin S 



L = 2x1*0 ~ cos 



(20) 

(21) 
(22) 



Let us now prove that the stability analysis can be 
reduced to two separate problems, 



L i ^ = A i V (i = 1,2), 



(23) 
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where the operators are 

Li = {L+ - L){L- + cos8), (24a) 
L 2 = [L- - cos S)(L+ + L) (24b) 

such that ImA 12 = and ReAi.2 > constitute neces- 
sary and sufficient conditions for the linear stability of 
the soliton. 

To this end, we notice that SU(4) rotation 








-1 





i 





1 







-1 





V- 


-1 





1 



provides for a unitary transformation, Cq 
with 



(25) 



pcp- 




Since the eigenvalues of £ and Cq coincide, we can 
consider the spectrum of the latter linear operator. Tak- 
ing into account that both C and Co are built of real 
coefficients, solutions can be looked for in the form of 
M 'i.2i u i2 ~ exp(iAz). Moreover, if A is an eigenvalue, 
then A is an eigenvalue as well (the overbar stands for 
the complex conjugate). In other words, the absence of 
an imaginary part of A, which is equivalent to the condi- 
tion that A = A 2 is real and positive, is a necessary and 
sufficient condition for the absence of the instability. 
As the next step, we consider the eigenvalue problem, 
Cq^> = A^f, where >J? = Pu, and we make use of the 



block structure of C\: 



- C 2 - 



( L 1 

4sin(<S)L_ L 2 

L 2 4sin((5)L + 

V Lx 



(27) 



where Lx l2 were introduced in Eqs. (|2~4"|) . Now, a straight- 
forward consideration demonstrates that A must coincide 
with either Ai or A2. Thus, the study of the stability of 
the dark solitons is reduced to eigenvalue problems (|23|) . 
Now, we notice that 



d 2 

-cos 8 = -^-2 - (xi +x){p 2 -u. 



'())■ 



La 



L 



L + 2(xi+x)uo- 



(28a) 
(28b) 



Therefore, taking into account Eq. (fTTj) . we conclude 
that the eigenvalue problem for operator L 2 is nothing 
but the standard stability problem for the black soliton 
in the defocusing medium, with the effective nonlinear- 
ity Xi + X- This problem is very well studied [TH, Il7j |. 
In particular, it is known that L + + L is positive defi- 
nite and L_ — cos (5 has only one negative eigenvalue and 
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FIG. 2: Regions of the stability and instability of the dark 
soliton (domains covered by dashed patterns and white ones, 
respectively) in the (b, xi) plane for fixed 7 = 0.7 and x = 1 
(a), and in the (6, 7)-plane for fixed xi — 2.5, x = 1 (b), or 
X = — 1 (c). In panel (a) inset presents the same regions as 
the main panel, but in the (p 2 ,Xi) pl ane iXi coincides with 
the left edge of the inset. 



one zero eigenvalue 16]. Moreover, it is known too [TtJ 
that the minimal eigenvalue of L 2 is positive. Thus, the 
eigenvalue problem for L 2 does not give instability, and 
our subsequent analysis is performed below for operator 
Lx, which may give rise to instability. 



V. NUMERICAL RESULTS 

The results of the numerical analysis of the linear sta- 
bility are depicted in Fig. [5] For the defocusing XPM, in 



the subdomain 



X < Xi 



< 



Xi (where Xi is the critical 



value, denoted in Fig. [2ja) by the blue vertical line), the 
dark-soliton stability region coincides with that for the 
CW background, which is — cos 5 < b < 2\i (cos 8) /(x — 
Xx) [or, equivalently, p 2 < (cos<5)/(x - Xi)] for Xi < X 
and b > — cos 8 for x < Xi < Xii see Figl2 a )- At the 
same time, in subdomain xi > Xi a dark-soliton's insta- 
bility "wedge" is present: as seen from Fig. the dark 
soliton is stable when — cos<5 < b < b\ or b 2 < b < oo. 
The value of the propagation constant, bx, at the lower 
edge of the "wedge" [the green line in Fig. HJa)] is almost 
independent of SPM coefficient xi (except for a small re- 
gion in a vicinity of the critical value Xi), while the upper 
edge, b = b 2 [the red line in Fig. Ufa)], is a quasi-linear 
function of xi- With the increase of 7 = sin 8 the dark- 
soliton instability "wedge" gradually shrinks [see Fig. 
[2jb)], disappearing at 7 = 1. For the focusing or zero 
XPM, with x = —1 or \ = 0, the dark soliton is stable 
at — cos 8 < b < bx , and unstable at b > bx , see Fig. [5Jc) . 
It is relevant to note that b\ does not depend on xi, and 
almost coincides with bx corresponding to the defocusing 
XPM [cf. Figs. d;b) and (He)], while b 2 coincides with 
the vertical line, xi = ~X- 

The linear stability analysis was completed by the di- 
rect simulations of Eqs. (JBJ). Typical examples of the per- 
turbed evolution of stable and unstable dark solitons are 
presented in Fig. [3] The predicted stability of the dark 
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FIG. 3: The evolution of field components |qi(:r,2:)| 2 (left 
column) and |g2(z,z)| 2 (right column) of dark soliton (|17|) 
with 7 = 0.7, xi = 2.5, x = 1, 6 = 0.5, p 2 w 0.347 (a), 
6 = 3.5, p 2 « 1.2 (b), or 6 = 8.5, p 2 » 2.633 (c). Panels 
(a), (b) and (c) correspond to points A, B, and C in Fig. [21 
respectively. 



soliton below the lower edge of the instability "wedge", 
i.e., at b < bi, is confirmed by the simulations. For 
the dark-soliton parameters corresponding to point A in 
Figs. [Ha) and [2(b), the evolution of fields component 
qi^(x,z) is shown in Fig. [3(a). Similarly, the evolu- 
tion of the dark soliton with parameters corresponding 
to point B, as well as the stability of the dark soliton 
with the parameters corresponding to point C in Figs. 
[3(a) and [2(b), are demonstrated in Figs. [3(b) and [3(c), 
respectively. 

For the focusing XPM, the stability and instability of 
the dark solitons [points A' and B' in Fig. [2(c)] is also 
confirmed by direct simulations, see Figs. 2(a) andH(b), 
respectively. 

The robustness of the dark solitons can be also be 
tested in interactions of two such solitons (kink-antikink 
pairs). Thus, in this case we use the initial condition at 
z = in the form of 



u(x) 



tanh )Jx±+x L- 1 - 




, (29) 



where t is the spatial separation between the two dark 
solitons. As can be seen from Fig. [SI in the VT- 




FIG. 4: The same as in FigO but for 7 = 0.5, xi = 2.5, 
X = -1, b = 0.5, p 2 « 0.910 (a), or b = 3.5, p 2 « 2.910 (b). 
Panels (a) and (b) correspond to points A' and B' in Fig. [2] 
respectively. 
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FIG. 5: The evolution of field components |gi(a;, z)\ (left 
column) and |q2(£, z)| 2 (right column) of dark-soliton pair 
pg)l with xi = 2.5, i = tt/2 and x = -1, 7 = 0.5, b = 0.5, 
p 2 w 0.910 (a), x = 1, 7 = 0.7, b = 0.5, p 2 w 0.347 (b), or 
X = 1, 7 = 0.7, 6 = 8.5, p 2 w 2.633 (c). Panels (a), (b) and 
(c) correspond to points A', A, and C in Fig. respectively. 



symmetric system the two dark soliton always [for both 
the focusing - Fig. [3(a) and defocusing - Figs. [5(b), [S(c) 
signs of the XPM] repel each other and start motion in 
opposite directions without self-destruction. The repul- 
sion from the boundaries of the x-domain in Fig. [5J hap- 
pens due to the implied periodic boundary conditions, 
and is equivalent to the repulsion between the dark soli- 
ton. As can be seen from the comparison of Figs. [S(b) 
and [5(c) , the increase of b (while separation I between 
the dark solitons is kept unchanged) results in reduction 
of the repulsion between the solitons, and, consequently, 
decrease of the solitons' "velocities" . The reason for this 
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FIG. 6: The evolution of field components |qi(:r,,z)| 2 (left 
column) and \q2(x, z)\ 2 (right column) of the dark soliton with 
initially separated components (|30[) . for D — tv/50. Other 
parameters are the same as in Fig. [5] 



phenomenon is that larger b corresponds to a smaller soli- 
ton width, and, as a result, a larger ratio of separation I 
to the soliton's width. 

Another possibility to set dark soliton in motion is to 
separate components in the initial condition, i.e., take 



(x, 0) = u(x + (-l)-' J D/2)e l (- 1 > J5 / 2 



(30) 



where u{x) is borrowed from Eq. (|17p . and D is the ini- 
tially imposed separation between the components. The 
results of the corresponding simulations are represented 
in Fig. [6] Here, for the defocusing XPM the "velocity" of 
the dark soliton does not strongly depend upon propaga- 



tion constant b [cf. Figs. 1307b) andHtc)]. The situation 
is completely different from the previous case [cf. Figs. 
\5{b) and[5jc)]. It should be noted that, in both Figs. 
[S] and [SJ the simulations were run for dark solitons with 
the propagation constants far enough from the stability 
margins 61.2 [points A', A and C in Fig. |2j. In the op- 
posite situation, for the dark solitons with propagation 
constants close to stability margins b±^, their motion may 
result in destruction, under certain conditions. 



VI. CONCLUSIONS 

To conclude we have reported the existence of stable 
vector solitons in the "PT-symmctric coupled nonlinear 
Schrodinger equation one of which has gain and another 
dissipation, whose strengths are equal. The found soli- 
tons have identical amplitude profiles but the phase dif- 
ference imposed by the gain-loss coefficients what ensures 
the balance between gain and losses. The stability of ei- 
ther backgrounds against which solitons propagate or of 
the solitons themselves are modified by dissipation and 
gain, what was confirmed by direct numerical simulations 
of the soliton propagation and interactions as well by the 
linear stability analysis. 
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